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We present a completed random matrix theory for staggered fermions which incorporates all taste 
symmetry breaking terms at their leading order from the staggered chiral Lagrangian. This is an 
extension of previous work which only included some of the taste breaking terms. We will also 
discuss the effects of taste symmetry breaking on the eigenvalues in the weak and strong taste 
breaking limits, and compare with some results from lattice simulations. 
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1. Introduction 

Staggered fermions are an inexpensive and popular way of discretizing quarks on a space-time 
lattice. They preserve a chiral symmetry with a local action at the expense of containing extra 
modes that quadruple the number of quark species being simulated. The quadrupled species, called 
tastes, are mixed at nonzero lattice spacing, but are expected to produce four independent quark 
flavors in the continuum limit. Evidence of this behavior can be seen in the low eigenvalues of 
the staggered Dirac matrix. At small lattice spacings the eigenvalues tend to cluster into distinct 
quartets which represent the four tastes. As the lattice spacing is decreased, the eigenvalues within 
a quartet will become more degenerate [[Tj]. 

In [||] a Random Matrix Theory (RMT) for staggered lattice fermions was introduced to de- 
scribe the low eigenvalues of the staggered Dirac operator. The staggered RMT (SRMT) adds addi- 
tional terms to the standard chiral Random Matrix Theory that have the appropriate symmetries of 
the lattice operator. These additional terms in the SRMT reproduce the known 0(a 2 ) (ignoring any 
extra factors of a s which we drop for convenience) taste breaking terms that appear in the staggered 
chiral Lagrangian. 

The SRMT is equivalent to the staggered chiral Lagrangian at zero momentum. This equiv- 
alence has been demonstrated for the fermionic partition function, but is only conjectured for the 
partially quenched case, which is related to the eigenvalues of the Dirac operator. Here we will di- 
rectly compare estimates of the size of the taste breaking in lattice simulations determined from the 
staggered chiral Lagrangian (in the /^-regime) with estimates obtained from comparing low eigen- 
values to the SRMT predictions (in the £-regime). Our initial tests show good agreement with the 
predictions of SRMT in support of the conjectured equivalence of the partially quenched theories. 

2. Staggered Chiral Lagrangian 

The effective chiral Lagrangian for staggered fermions at order a 2 is given by [||, [|] 



where (X) stands for the trace of X, and F and Zo are the low energy constants (LECs) related to 
the pion decay constant (with the convention that the physical value for F rj 131 MeV) and the 
magnitude of the chiral condensate, respectively. The taste breaking terms can be divided into two 
parts ~f = ~f\ t + % f . The first part contains the single-trace terms (with ^ = y^) 




(2.1) 




(2.2) 



and the second part has the two-trace terms 



- % t = C 2 v\X[{^U) (^U) + h.c] + C U W [(feU) (^U)+h.c] 




(2.3) 
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3. Staggered Chiral Random Matrix Theory 

The (fermionic) staggered chiral random matrix theory partition function can be written as 

Ni 



Z S rmt= dWpo{W)p T (T)Yldct{D + m f ) (3.1) 

J /=i 



with 



/ iW \ 

D= [iw< )® l4+aT - M 

where W is a (N + v) x N complex matrix with v the absolute value of the topological charge and 
T incorporates the taste breaking terms. The matrix potential for W is conveniently a Gaussian, 

Po (W) = exp(-aN(W f W)) (3.3) 

with yfa = LqV/2N (V is the four volume). 

The taste breaking contribution to the SRMT (T) has eight terms that correspond directly with 
the eight taste breaking terms of the chiral Lagrangian. Its complete form was given in [||]. As an 
example, the Ca, term corresponds to 

I a n \ 

(3.4) 

where and are Hermitian matrices of size (N + v) x (N + v) and N xN, respectively. For 
convenience we can choose a Gaussian weight function for these matrices, 





One can then show that the (fermionic) RMT with this extra matrix term is equivalent to the zero- 
momentum staggered chiral Lagrangian with just a C4 taste breaking term. 

Note that while the correction to the RMT enters at order a, this still reproduces a term of order 
a 2 in the chiral Lagrangian. This is due to the taste breaking terms in the SRMT being traceless. 
As demonstrated in [||], when expanding the determinant of the SRMT Dirac matrix, the 0(a) term 
vanishes for this reason, and results in a partition function that has a leading corrections at 0(a 2 ) 
even though the SRMT Dirac operator has terms of 0(a). 

The two-trace terms can be incorporated in the SRMT in a couple of ways. One is to add terms 
such as 



±b tl I N 



®Zn5 (3.6) 



with a Gaussian weight for the scalar b^. This will give a contribution to the Cia and Csa terms. 
While this term will reproduce the correct term in the chiral Lagrangian, it has no analogue in the 
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lattice Dirac matrix. This term is of the form of a fluctuating taste-dependent mass, which is not 
present on the lattice. 

An alternative form for the two-trace terms is to simply modify the potentials for the matrix 
terms corresponding to the one-trace terms. By replacing the simple Gaussian weight with one that 
also includes (A^ ) , (B^ ) and (A^ ) (B^ ) in the exponential, the coefficients can be tuned to give 
the correct two-trace terms in the chiral Lagrangian [§]. 

4. Generalized Staggered Random Matrix Theory 

If we look at the final form of the SRMT Dirac matrix, we see that it is the most general ma- 
trix that is consistent with the symmetries of the staggered Dirac operator; it is anti-Hermitian and 
anticommutes with the staggered chiral symmetry matrix 75 (g) ^5. As mentioned in the previous 
section, the single-trace terms from the chiral Lagrangian are reproduced by considering indepen- 
dent Gaussian weights for the remaining matrix elements. Meanwhile the two-trace terms can be 
reproduced by adding two-trace terms to the RMT weights. 

In this manner one could consider generalizing the SRMT to include more terms in the weight 
function in an attempt to reproduce higher order terms in the chiral Lagrangian (at zero momen- 
tum). One can then draw an analogy between the formulation of a RMT and of an effective La- 
grangian. For the Lagrangian one includes all terms, up to some order, that are consistent with the 
symmetries. Likewise for the RMT one can consider a matrix containing all elements consistent 
with the symmetries of the Dirac operator, with a generalized weight function for these matrix ele- 
ments up to some level of complexity. Of course the mapping from the RMT potential to the chiral 
Lagrangian at higher order may not be as simple as in the SRMT presented here, but one could 
speculate that a generalized RMT could reproduce all higher order terms of the chiral Lagrangian. 
Again this equivalence only holds for the zero momentum Lagrangian, but it is also possible to 
reduce the full Lagrangian to an effective zero momentum Lagrangian (for a recent example see 
[§]), which then might be representable as a RMT. 

5. Dominant form of taste breaking 

Most of the taste breaking coefficients can be measured in lattice simulations by comparing to 
staggered chiral perturbation theory (S#PT) results. The four single-trace parameters are uniquely 
determined from the splittings of the pion masses at leading order [Q]. Typically the C4 term is 
found to be the dominant contribution to the pion spectrum [||]. The two-trace terms enter in 
S^PT formulas in the combinations Cy A = (C 2 {v,a} =tC 5 {v,A})/2- They don't contribute to the pion 
splittings at leading order, but the C A and Cy terms do appear in one-loop expressions, while the 
Cj~ and Cy terms do not [Q]. In lattice simulations C A has been found to be larger than Cy [§]. 

Both of the dominant terms, C4 and C A , come from same term in the SRMT with the more 
general form for the weight function. This supports the idea of constructing the SRMT from a 
single matrix with the correct symmetries and with a generalized weight function. Based on the 
lattice measurements, the leading contribution to taste breaking in the staggered Dirac matrix has 
an axial- vector taste structure. If one imagined rotating the staggered Dirac matrix into a taste basis 
and then expanding in powers of a, the dominant correction at order a would then have the same 



4 



Chiral random matrix theory for staggered fe unions 



James C. Osborn 




form as in (^4). Of course the exact potential for the lattice Dirac matrix would be much more 
complicated than in the SRMT, but the leading effects at low energy can be captured by the SRMT 
weight function considered here. Among the terms in the SRMT potential corresponding to the 
axial-vector matrix, we have no reason to favor one over the other, so we would naively expect 
them to be of similar magnitude. We would then expect the corresponding terms in the chiral 
Lagrangian to be of similar order, and also dominant over the corresponding terms with other taste 
symmetries, which is consistent with lattice measurements. 



6. Extracting LECs from RMT 

The RMT predictions for the eigenvalue correlations can be used to extract LECs from lattice 
simulations. For example Lq can be obtained by fitting to the eigenvalue density. Additionally F 
can be obtained from the correlations of eigenvalues with an imaginary chemical potential [0]. 

If the taste breaking is small enough then these methods can apply directly to staggered eigen- 
values by replacing each quartet of eigenvalues with its average [§]. In this case one could also try 
to extract the taste breaking parameters from the splittings of the eigenvalues within the quartet. In 
practice, it would likely be too difficult to extract all the parameters, but if one assumes that the C4 
term is dominant, then it should be possible to estimate it from comparison to the SRMT. 

However, if the taste breaking is large, then the higher order taste breaking terms can become 
important. In this case the effective chiral Lagrangian reduces to a single flavor for the remaining 
staggered chiral symmetry with a new set of LECs that are in principle unrelated to the original 
ones [Q]. Thus extracting the LECs from low eigenvalues when there aren't clear quartets present 
may not yield the continuum LECs in chiral Lagrangian. 
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Figure 2: Individual integrated densities of the lowest four eigenvalues from the 32 4 lattice ensemble and 
the SRMT. The dominant taste breaking parameter in the SRMT is set to a 2 VC 4 = 0.3 (left) and 0.2 (right). 

7. Comparison to lattice simulations 

As an initial test of the SRMT, we will compare the leading order taste breaking term obtained 
from fitting the Dirac operator eigenvalues to the SRMT with that obtained from fitting the pion 
spectrum to the staggered chiral Lagrangian. Since it is difficult to get a single lattice ensemble 
that we could use for both measurements we will use two ensembles with all parameters identical 
except for the volume. 

For the pion masses we use an ensemble from the MILC collaboration 2+1+1 flavor HISQ runs 
[@]. The ensemble has a volume of 32 3 x 96 with a lattice spacing of a « 0.09 fm and with a light 
quark mass m/ = m s /5. From the pion mass splittings we can get the single-trace taste breaking 
terms in the chiral Lagrangian. Using a value of F = 131 MeV, this gives 

a 2 VC x = 0.03(8), a 2 VC 3 = -0.03(4), a 2 VC 4 = 0.84(4), a 2 VC 6 = 0.03(3) . (7.1) 

We can see here that C\ is clearly dominant, as expected, and that the other coefficients are all 
consistent with zero. 

We generated a new ensemble of 430 lattices of size 32 4 with all other parameters the same as 
the previous ensemble. From the volume scaling we expect to find 



a 2 VC 4 



0.28(1) 



(7.2) 



on this new ensemble. We then compare lattice eigenvalues with numerical simulations of the 
SRMT with only the C4 taste breaking term. For this comparison we choose to use the number 
variance statistic. This shows the fluctuations (variance) in the number of eigenvalues in an inter- 
val starting at zero versus the average number in the interval. We evaluated it numerically from 
simulations of the SRMT with N = 400. 

In Figure[I] we P^ ot tne num ber variance of the lattice eigenvalues against the SRMT at different 
values of a VC4. We can see that the best agreement for small intervals (up to around 2 to 4 



eigenvalues) is at a 2 VC 4 0.3. This is in good agreement with our prediction (7.2) obtained from 
the pion mass splittings. As the length of the interval grows, the lattice results start to move away 
from the SRMT result. This is likely due to higher momentum modes entering on the lattice, that 
aren't captured in the RMT. This happens at the QCD equivalent of the Thouless energy [Q], which 
for this ensemble, in units of the average eigenvalue spacing, is F 2 \fV « 3.5. This is consistent 
with our observations from the number variance. 
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As a further check that the SRMT describes the low eigenvalues of the staggered Dirac oper- 
ator, we look at the individual integrated densities of the lowest four eigenvalues. In Figure || we 
see the integrated density from the 32 4 lattice along with that from the numerical simulations of the 
SRMT at two different values of the taste breaking parameter a 2 VC^ = 0.2,0.3. We can see that 
the value of 0.3 fits the lattice data much better than at 0.2, again confirming that the predictions 
of the SRMT fit the lattice data well and give estimates of the taste breaking parameters that are 
consistent with the staggered chiral Lagrangian. 

8. Summary 

We have shown a chiral RMT that incorporates all leading order taste breaking terms from 
staggered chiral Lagrangian. The SRMT can be constructed by considering a RMT with same 
symmetries as the staggered Dirac operator and a generalized weight function. Initial tests show 
that the predictions of the SRMT are in good agreement with lattice simulations within the range 
of validity of the SRMT. Additionally, using the predictions of the SRMT, the dominant taste 
breaking parameter can be extracted from the low eigenvalues of the staggered Dirac operator and 
gives consistent results with that obtained from the pion mass splittings. 
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